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ABSTRACT 

Anomalous microwave emission (AME) has been observed in numerous sky regions, 
using different experiments in the frequency range ~ 10 — 60 GHz. One of the most 
scrutinized regions is G159.6-18.5, located within the Perseus molecular complex. In 
this paper we present further observations of this region (194 hours in total over « 
250 deg^), both in intensity and in polarization. They span four independent frequency 
channels between 10 and 20 GHz, and were gathered with QUIJOTE, a new GMB 
experiment with the goal of measuring the polarization of the GMB and Galactic 
foregrounds. When combined with other publicly-available intensity data, we achieve 
the most precise spectrum of the AME measured to date in an individual region, 
with 13 independent data points between 10 and 50 GHz being dominated by this 
emission. The four QUIJOTE data points provide the first independent confirmation 
of the downturn of the AME spectrum at low frequencies, initially unveiled by the 
COSMOSOMAS experiment in this region. We accomplish an accurate fit of these 
data using models based on electric dipole emission from spinning dust grains, and 
also fit some of the parameters on which these models depend. 

We also present polarization maps with an angular resolution of « 1° and a 
sensitivity of « 25 /rK/beam. From these maps, which are consistent with zero po¬ 
larization, we obtain upper limits of H < 6.3% and < 2.8% (95% C.L.) respectively 
at 12 and 18 GHz, a frequency range where no AME polarization observations have 
been reported to date. These constraints are compatible with theoretical predictions 
of the polarization fraction from electric dipole emission originating from spinning 
dust grains. At the same time, they rule out several models based on magnetic dipole 
emission from dust grains ordered in a single magnetic domain, which typically pre¬ 
dict higher polarization levels. Future QUIJOTE data in this region may allow more 
© 2014 RAS stringent constraints on the polarization level of the AME. 
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1 INTRODUCTION 

Now that the observations of the temperature anisotropies 
of the Cosmic Microwave Background (CMB) have reached 
levels of sensitivity and angular resolution that have allowed 
the determination of the main cosmological parameters with 
accuracies close to 1% ( Bennett et al.|2013 Planck Collabo- 
ration et al. |2014c[ ), attention has shifted to the study of the 


polarization of these anisotropies. CMB polarization also en¬ 
codes a wealth of cosmological information. In particular, it 
is of paramount importance to study the B-mode anisotropy, 
which can only be created by vector or tensor perturba¬ 
tions, such as those due to primordial gravitational waves 
I Kamionkowski et al.|[T9^ Zaldarriaga fc Seljak||19^ l. As 
all inflationary scenarios are predicted to give rise to gravi¬ 
tational waves, this B-mode signal is a definite test for infla¬ 
tion. Following previous upper limits placed by experiments 


like QUIET (QUIET Collaboration et al. 20121 or BICEP 


(BICEPl Collaboration et al. 20141, the first detection of 


this signal was claimed in data from the BICEP2 experi¬ 
ment at 150 G Hz, with a level of the tensor-to-scal ar ratio 
r = 0.20i'lg Q 5 (BICEP2 Collaboration et al. 20141. How¬ 
ever, soon after these results were published, some papers 


(Flauger et al. 2014 Mortonson & Seljak 20141 appeared 


pointing out that the BICEP2 team could have underesti¬ 
mated the contribution from polarized thermal dust emis¬ 
sion, leading to a wrong interpretation of their signal. This 
has been recently supported by data from the Planck satel¬ 
lite, which have shown that the level of polarized dust emis¬ 
sion in the region of the sky covered by BICEP2 could form 


a significant component of the measured signal ( Planck Col 
[laboration et al!]|2014b[ ), thus causing a likely reduction in 
the level of cosmological signal that can be inferred from the 
BICEP2 results. 

It is therefore clear that, any unambiguous B-mode de¬ 
tection requires a full assessment of the level of foreground 
contamination and ideally confirmation by independent ex¬ 
periments operating at different frequencies. QUIJOTE is 
one of these experiments, which operates in the frequency 
range 10-40 GHz, and is subject to different systematics and 
foregrounds. QUIJOTE data will also be useful for measur¬ 
ing the polarization of low-frequency foregrounds, which are 
known to dominate the underlying primordial B-mode po¬ 
larization over a wide frequency range, and therefore need 
to be characterized and removed from the observed maps. 
Synchrotron emission from relativistic electrons is the most 
important contaminant at low frequencies in polarization. 
In intensity, free-free emission and the so-called anomalous 
microwave emission (AME) also show up in this frequency 
range. While the former is known to be practically unpolar¬ 
ized ( |Rybicki &: Lightman|jl97^ , very little is known about 
the polarization level of the AME. Since its discovery in the 
90s I Kogut et al.[1996a|b Leitch et al.||1997 1, many observa¬ 


tions in large-sky areas Oliveira-Costa et al.||199^|1999[ 

Davies et al. [200^1 and in individual Galactic (Finkbeiner et 


al.|2002| [Watson et al.|2005||Casassus et al.|2006[ Dickinson 

eTakf^TOI XMr'Con^rtium'er'ai‘1200^ Ubb^etaLl^oT^ 


Genova-Santos et al.|20lT] [Vidal et al.|2011||Planck Coll^ 


oration et al.|2011 2014aI and extragalactic (Murphy et al. 


20101 clouds have contributed to the understanding of the 
physical properties of this emission. A great deal of effort has 
also been dedicated to theoretical modelling of AME. Elec¬ 
tric dipole radiation from very small and rapidly rotating 


dust grains in the interstellar medium (Draine & Lazarian 

1998[ [Ali-Hai'moud et al.[[2009 

[Hoang et al.[[2010| [Ysard & 

Verstraete[[2010 Silsbee et al. 

2011 

1 , the so-called spinning 


dust emission, is the scenario that best fits the observations. 
An alternative mechanism based on magnetic dipole emis¬ 
sion has also been proposed (Draine fc Lazarian|1999 1, with 


a spectrum peaking at higher frequencies. 

The characterization of AME polarization properties is 
important not only to get further insight into its physical 
mechanism, but to assess how it could affect experiments 
pursuing the primordial B-mode signal in the frequency 
range 10 to 50 GHz. There are some theoretical studies of 
the AME polarization in the literature. After the first pre¬ 


dictions by Lazarian & Draine (20001 of spinning dust polar¬ 


ization, more recently Hoang et al. (20131 presented a model 


based on observations of the UV polarization bump, accord¬ 
ing to which the maximum polarization fraction would be 
« 3%, at a frequency of 5 GHz. In regard to the magnetic 
dipole emission model, while Draine fc Lazari^ ( [1999 1 pre¬ 
dicted high polarization fractions (up to 40 %) in the case 
of dust grains with atomic magnetic moments oriented in a 
single domain, Draine & Hensley ( |2013| recently presented 
a more realistic model, with randomly-oriented magnetic in¬ 
clusions, which results in lower polarization degrees (< 5% 
in the range 10—20 GHz). However, not much is known, from 
the observational standpoint, about the polarization proper¬ 
ties of the AME. Using WMAP 5-year data [Macellari et al.[ 
(20111 set an upper limit H < 5% (this and other upper lim¬ 
its that will be referred to in this section are at the 95% confi¬ 
dence level) on the polarization fraction of the diffuse AME. 
Some other constraints refer to individual clouds. [BattistelE] 


et al. (20061 observed the Perseus molecular complex with 
the COSMOSOMAS experiment and derived H = 3.4]_i.9% 


at 11 GHz. Casassus et al. (20081 reported an upper limit 


of < 1.0% at 31 GHz on the p Ophiuchi molecular cloud 
using the Cosmic Background Imager, whereas [Mason et al.[ 
(20091 found a maximum of < 2.7% with the Green Bank 


Telescope at 9 GHz. More recently, [Lopez-Garaballo et al.[ 
| |2011[ ) obtained an upper limit of < 1.0% at 23 GHz on 
the Perseus molecular complex using WMAP 7-year data. 


Shortly after this paper Dickinson et al. (20111, using the 


same data, obtained < 1.4% in the same region, and < 1.7% 
in p Ophiuchi. A detailed review of all these observations, 
plus some updated constraints in some regions, have been 
presented in [Rubino-Martm et al.| ( |2012b| . 

In this paper we present the first results obtained with 
the QUIJOTE experiment, that are based on observations 
of G159.6-18.5 in the Perseus molecular complex, one of the 


most studied AME regions in the sky ( 

Watson et al. 

2005 

Tibbs et al.[[2010[ [2013[ [Planck Gollaboration et al. 

2011 

2014a[). QUIJOTE observations cover the frequency range 
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10 — 20 GHz, where only the COSMOSOMAS experiment 
had provided observations of the AME before. The goal of 
this paper is twofold: to confirm the downturn of the AME 
spectrum at frequencies below 23 GHz through similar spec¬ 
tral sampling but completely independent results to those 
provided by GOSMOSOMAS, and to set constraints on the 
polarization level of the AME in the so far unexplored spec¬ 
tral region between 12 and 20 GHz. Section 2 is dedicated 
to the description of the observations and the basics of the 
data reduction. Our main results are presented in section 3, 
while the conclusions of this work are discussed in section 4. 


2 DATA AND METHODOLOGY 

2.1 QUIJOTE data 

The new data presented in this article were acquired with the 
QUIJOTE experiment. QUIJOTE is a collaborative project 
that consists of two telescopes and three polarimeter instru¬ 
ments covering respectively the frequencies 10-20, 30 and 
40 GHz, and located at the Teide Observatory (2400 m a.s.l.) 
in Tenerife (Spain). The main science driver of this exper¬ 
iment is to perform observations of the GMB polarization 
to constrain the B-mode signal down to r = 0.05, and to 
characterize the polarization of low-frequency foregrounds, 
mainly synchrotron emission and the AME, so that this sig¬ 
nal can be removed from the primordial maps to a level that 
will permit reaching the previous level of r. The two QUI¬ 
JOTE telescopes are based on an offset crossed-Dragone de¬ 
sign, with projected apertures of 2.25 m and 1.89 m for the 
primary and secondary mirrors, and provide highly sym¬ 
metric beams (ellipticity > 0.98) with very low sidelobes 
(?; —40 dB) and polarization leakage (^ —25 dB). The first 
instrument to be fielded on the QUIJOTE first telescope is 
a multi-frequency instrument (MFI) with 4 horns covering 
the frequency range 10 — 20 GHz, and with angular reso¬ 
lutions close to one degree. These detectors are fitted with 
MMIG low-noise amplifiers (noise temperature better than 
10 K), and use stepped polar modulators to measure the po¬ 
larization of the incoming radiation, providing instantaneous 
sensitivities of « 650 fj,K in four individual frequencies 
with nominal figures: 11, 13, 17 and 19 GHz. The median 
integrated PWV above the observatory is « 4 mm, giving 
a zenith atmosphere temperature of ~ 2 K at 11 GHz and 
~ 5 K at 19 GHz. The MFI saw first light on November 
2012 and since then it is performing routine observations of 
different Galactic and cosmological regions. The second in¬ 
strument consists of 31 polarimeters at 30 GHz (TGI, thirty- 
gigahertz instrument), and is based on the same design of 
the MFI except that the polarization modulation is achieved 
electronically through phase switches. This instrument will 
be commissioned during 2015. Finally, the third instrument 
is planned to have 31 polarimeters at 40 GHz (FGI, forty- 
gigahertz instrument). Using the TGI and the FGI, which 
will provide instantaneous sensitivities of 50 s^^^, we 

plan to survey an area of 3000 deg^ down to a projected 


sensitivity of 1 /rK/beam. A more detailed description of 
the technical and scientific aspects of this project can be 
found in Rubino-Martfn et al. (2012a I or in Rebolo et al. 
(in preparation). 


2.1.1 Observations 

The observations covering the Perseus molecular complex 
were carried out between December 2012 and April 2013 us¬ 
ing the MFI, in four frequency bands centred at 11.2, 12.9, 
16.7 and 18.7 GHz. The beam FWHMs are 0.87°, for the 
two lower frequencies, and 0.65°, for the two high frequency 
bands. The observations consisted of raster scans at constant 
elevation in order to minimize the effect of atmospheric vari¬ 
ations. Each scan had an amplitude in azimuth direction of 
12 degrees centred around the coordinates RA = 3’’52"', 
Dec. = 34°. This position was chosen to be equidistant be¬ 
tween the AME cloud G159.6-18.5 and the California HII 
region (NGC1499), which is also observed in the scans and 
is used as a null test for zero polarization. In each raster 
scan, of total integration time of 30-35 minutes, the telescope 
moves back and forth in azimuth at a velocity of 1 deg/s. A 
total of 336 raster scans were performed, in four positions 
of the polar modulators (0°, 22.5°, 45° and 67.5°) in order 
to minimize systematics. In total 194.4 hours of data were 
accumulated, 23% of which were removed due to being af¬ 
fected by bad weather or instrumental effects, resulting in 
a final effective observing time of 148.9 hours. As the MFI 
horns are separated typically by 5 degrees on the sky, the 
sky coverage of each one is different. The total sky area cov¬ 
ered by the four horns was respectively 176, 184, 277 and 
261 deg^. 


2.1.2 Amplitude calibration 


We determine the gain calibration factors (giving the conver¬ 
sion from voltage measured in the detectors to temperature 
on the sky) using intensity measurements on Gas A, which 
is ideal as it is bright and has a very low degree of polar¬ 
ization. We do daily 25-min raster scans of 10 x 10 degrees 
around this source from which we derive the gain calibration 
factors for each channel. The output signal of each channel 
is a combination of the three Stokes parameters I, Q and 
U. According to WMAP the polarized intensity of Gas A 
at 22.8 GHz is P = 0.81 ± 0.05 Jy (polarization fraction 
H = 0.35±0.02%) (Weiland et al. 20111, which is low enough 
not to be detected in a single raster scan. We can therefore 
safely assume Q = U = 0, and use the I flux to calibrate. 
To achieve this we use the modelled spectrum of | Weiland| 
et al. (20111, which is obtained by fitting a combination of 


WMAP 7-year data and other ancillary data to a logarith¬ 
mic quadratic function. This function is then integrated over 
the measured spectral transmission of each frequency band 
to obtain the reference Gas A fluxes associated with each 
channel. 

Finally, as the previous model is referred to epoch 2000, 
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in order to account for the secular decrease of the Cas A 
flux (typically 0.5%/year), we use the Hafez et al. (20081 
model, which was derived using VSA observations, in order 
to refer the final fluxes to the time of the observations. In or¬ 
der to circumvent possible uncertainties associated with this 
secular variation, a more suitable calibration source would 
be Jupiter. However, owing to its small angular size, this 
source is severely diluted in the QUIJOTE beams and a 
large number of observations would be required for it to be 
used as primary calibrator. We point out that, in any case, 
the gain calibration will only affect the modelling of the SED 
in G159.6-18.5. The gain calibration factors will cancel out 
when dividing the polarized intensity by the total intensity, 
and therefore the inferred polarization fractions, which are 
one of the main goals of this paper, are insensitive to the 
absolute flux calibration. 


2.1.3 Polarization calibration 


One of the main steps of the data processing is the calibra¬ 
tion of the polarization angle ipo, which is defined as the 
reference position angle of each polar modulator. To accom¬ 
plish this we use Tau A (also known as the Crab Nebula) as 
a calibrator, which is the brightest polarized source in the 
sky in the microwave range. We perform daily 25-min raster 
scans of 10 x 10 degrees around this source, from which we 
derive a polarized flux that is a function of the intrinsic Q/1 
and U/I, of the position of the modulator relative to ipo and 


of the parallactic angle (pp. WMAP 7-year results (Weiland 
et al.|20TT]) show that the Q/I and U/I ratios for Tau A do 


not significantly vary (less than 2%, which is consistent with 
the error of the measurement) between 23 and 94 GHz. We 
therefore assume that these factors will remain equally un¬ 
changed down to 10 GHz, and use as reference the WMAP 
measurements at 22.8 GHz, the closest frequency. As we also 
know ifip, we can therefore fit for (fio- 

Using 191 raster scans on Tau A throughout a year we 
have checked that the recovered polarization angle, ipo, is 
stable over time. We then combine all these observations to 
derive a unique value for each horn. The accuracy on the de¬ 
termination of this angle is respectively 0.8° and 0.4° for the 
two horns that will be used in the polarization analyses that 
will be presented in this paper. In QUIET, an experiment 
which has similarities with QUIJOTE, a precision of 0.5° 
is achieved by using a combination of Tau A observations 


with a sparse-wire-grid calibrator (QUIET Collaboration et 
|al.||2012| . We point out however that the accurate determi¬ 
nation of this angle is important only to derive precise Q and 
U fluxes. An incorrect angle will result in a mixing of flux be¬ 
tween Q and U, but the polarized intensity P — 
will remain unchanged. As in this analysis we will get con¬ 
straints on the polarization fraction, the accurate calibration 
of the polarization angle is unimportant. 


2 . 1.4 Map making 

The four output channels of each frequency band contain a 
combination of three Stokes parameters I, Q, U. The sum 
of these channels, after calibration of their individual gains, 
give I, while the subtraction of pairs of channels gives the 
following combination of Q and U: 

Uaub = Qsin(4yjpm + 2(/Pp)-I-f/cos(4(ppm-I-2v5p) , (1) 


where ippm is the position angle of the polar modulator, 
whose reference position is calibrated following the proce¬ 
dure explained in section |2.1.3[ and ipp is the parallactic 
angle. Out of the four channels, two are correlated and there¬ 
fore are affected by the same 1// noise, while the two other 
are not. To reconstruct the polarization signal we then only 
use the two correlated channels, in order to minimize the 1// 
contribution. The typical knee frequencies of our receivers 
are /k ~ 10 — 40 Hz depending on the channel. However, for 
the measurement of polarization, the subtraction results in 
much lower values of /k ~ 0.1 — 0.2 Hz. In order to further 
reduce the 1/f noise in the final maps, we apply a filter on 
the time-ordered-data (TOD) by subtracting the median of 
the data in intervals of 20 s, after binning the data at 50 ms. 

Under the assumption that the filtered TODs are dom¬ 
inated by white noise, we consider the noise covariance ma¬ 
trix to be diagonal, a hypothesis that considerably simplifies 
the map-making. As the response of our instrument to po¬ 
larization is a combination of Q and U, in order to recover 
these Stokes parameters in each pixel we have to combine 
all the samples lying in that pixel corresponding to different 
angles g} = + 2(^p. To do so we use two independent 

strategies. The first one consists on producing 100 maps, 
each one corresponding to g) angles within a given bin, and 
then using the 100 values of each pixel to find the best-fit 
solution for Q and U from equation The second strategy 
builds on an analytical minimization. The different pa¬ 
rameters, which are combinations of sines and cosines of the 
g} angles, are grouped in each pixel, and at the end of the 
process the Q and U values are computed using the ana¬ 
lytical formulas that result from this minimization. In both 
cases the data samples are weighted according to their noise, 
which is calculated from the standard deviation calculated 
during the binning of the TODs. 

To produce the final maps we use a HEALPix pixeliza- 
tion ( Gorski et al.|2005 l with A^side = 512 (pixel size 6.9 ar- 
cmin), which is sufficient given the beam FWHM. While we 
have checked that the maps resulting from the two strate¬ 
gies are almost identical, in the subsequent analyses we use 
those resulting from the second method. 


2.2 Ancillary data 

All the polarization data that will be used in this article 
comes from the QUIJOTE experiment. However, in order to 
obtain the full spectral energy distribution (SED) of G159.6- 
18.5, from which the residual AME fluxes will be inferred, 
we use ancillary data from other experiments. In the low- 
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frequency range we use the 

Haslam et al. 

(1982 

1 maprl at where ni is the number of pixels in the aperture, and T and 

0.408 GHz, the 

Berkhuijsen 

(19721 ma|;[^[at 0.820 GHz and Tj represent respectively the pixel thermodynamic temper- 

the Reich & Reich (19861 map at 1.4 GHz. 

atures in the aperture and in the background annulus. The 


At 10.9, 12.7, 14.7 and 16.3 GHz, similar frequencies 
to QUIJOTE, we use data from the COSMOSOMAS ex¬ 


periment (Watson et al. 20051. In order to minimize the 


1// noise, the data from this experiment was filtered by the 
suppression of the first seven harmonics in the FFT of the 
circular scans, which results in a flux loss on large angular 
scales. For this reason, the comparison with other experi¬ 
ments that preserve all the angular scales is not straight¬ 
forward. It is necessary to account for the flux lost, as it 


was done in Planck Collaboration et al. (20111. The fluxes 


presented in that paper are already corrected, so we directly 
take those fluxes. 

We al so use data from the 9-year release of the WMAP 
satellit^ ( [Bennett et al.|[2013[ ) to provide flux estimates at 
frequencies 23, 33, 41, 61 and 94 G Hz. Recent data from the 
first release of the Planck missiorj^ ( jPlanck Collaboration et' 
al. 2014c I cover the frequencies 28, 44, 70, 100, 143, 217, 


353, 545 and 857 GHz. We also download the released Type 
1 CO maps ( Planck Collaboration et al.||2014e l, which are 
used to correct the 100, 217 and 353 GHz frequency maps 
from the contamination introduced by the CO rotational 
transition lines (1-0), (2-1) and (3-2), respectively. Finally, 
in the far-infrared spectral range we use Zodi-Subtracted 


Mission Average (ZSMA) COBE-DIRBE maps (Hauser et 
^[T9M| at 240 /rm (1249 GHz), 140 /rm (2141 GHz) and 
100 (2998 GHz). 


2.3 Methodology for flux estimation 

Intensity and polarization fluxes will be calculated by apply¬ 
ing an aperture photometry integration on the maps. This 
is a well known and widely used technique in this context 
(jLopez-Caraballo et al.|2011| [Dickinson et al. [20111 |Genova-| 


Santos et al. 20111, consisting of integrating temperatures 


of all pixels within a given aperture, and subtracting a back¬ 
ground level calculated in an external ring. Instead of the 
mean, following Planck Collaboration et al. (20111, we chose 


to use the median of all the pixels in the external ring as the 
estimate of the background level. The median is a better 
proxy for the real level in cases of strongly variable back¬ 
grounds with many outlier pixels. The flux is then given by 


Si, = a{v) 


E nl rji 
1=1 


-Ti 


m 


( 2 ) 


^ We use the map supplied by 
^ We projected the 


(20031 


t 


Platania et al. 

map downloaded from 
http://www.mpifr-bonn.mpg.de/survey.htinl into HEALPix 
pixelization. 

^ Downloaded from the LAMBDA database, 
http://lambda.gsf c.nasa.gov/ 

^ Downloaded from the Planck Legacy Archive, 
http: //www. sciops. esa. int/index. php?pro j ect=plcuick.i£page= 
Planck_Legacy_Archive. 


median is calculated over the 712 pixels in this annulus. The 
function a{i') gives the conversion factor from temperature 
to flux, 


a{u) = 




sinh 


-2 / hv \ 
V 2fcbT'cmb / 


niHp 


( 3 ) 


where h and kh are the Planck and Boltzmann constants, 
c the speed of light, Tcmb = 2.725 K (Mather et al. 19991 
the CMB temperature and Hpix the solid angle subtended 
by each pixel. 

The determination of the error associated with the pre¬ 
vious estimate is crucial for the results of this paper. In 
a hypothetical case of perfect white uncorrelated noise, it 
could easily be estimated through: 


a-{SA = aiy) 


E nl 2 2 

i = l I 

n? 2 n‘i 


1/2 


( 4 ) 


where a represents the error of each pixel. 

However, in QUIJOTE the instrument noise is corre¬ 
lated due to the presence of 1// residuals, and also the galac¬ 
tic background fluctuations introduce, mainly in intensity, 
an important contribution to the error which is correlated 
on the order of the beam size. Ideally we should then use 
the covariance matrices of the instrument and background 
noises. The former can be extracted through a characteriza¬ 
tion of the 1// noise spectrum, however the latter is difhcult 
to determine. Instead, in the previous equation we can intro¬ 
duce in the denominator the number of independent pixels in 
the aperture and in the ring, which we will denote respec¬ 
tively as n'l and n^. The pixel variance will be calculated 
from the pixel-to-pixel standard deviation of all the pixel 
temperatures Tj in the background, a(Tj). Obviously, the 
standard deviation of the pixels in the aperture would be 
biased by the presence of the source. On the contrary, the 
standard deviation of the pixels in the ring gives a reason¬ 
able estimate of the contributions of the background and of 
the instrumental noise to the true error. Therefore, the final 
equation that we will use to estimate errors in this paper 
reads as: 


u(S.) = a(^)a(T,)J^ + |^ 


( 5 ) 


In the case of the error being completely dominated by 
the background, then we could use for n( and n '2 the num¬ 
ber of beams in the aperture and in the background, rii 
and n 2 . However, while being particularly strong in inten¬ 
sity, the background fluctuations from the Galactic emission 
are not so important in polarization. For this reason, in this 
case the relative contribution from the 1// residuals to the 
uncertainty is significant. To quantify this, we selected 20 
random positions around our source, G159.6-18.5, and per¬ 
formed flux integration on those positions using the same 
aperture and ring sizes. The standard deviation of those 
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Figure 1. QUIJOTE intensity map at 11 GHz of the whole region covered by the observations (left), in comparison with the WMAP 
9-year map at 23 GHz (right). The positions of the G159.6-18.5 molecular cloud, the California HII region (NGC1499) and the 3C84 
quasar are marked with crosses. The QUIJOTE map encompasses 277 deg^, contains in total 149 h of observations, and its RMS is 
80 /rK/beam. By comparing the relative amplitudes of California and G159.6-18.5 it can easily be noted how the presence of AME 
boosts the brightness of G159.6-18.5 at 23 GHz. 


fluxes gives a reasonable estimation of the true noise of our 
flux estimate. From this analysis we determined that for in¬ 
tensity n \^2 = while for polarization n'^ 2 = 2ni 2 . This 
is what we will use in our estimation of the flux errors. 


In cases of low signal-to-noise fluxes, or when placing 
upper limits on the polarized flux P, as it will be our case, 
it is necessary to de-bias the fluxes derived from the aper¬ 
ture photometry integration. This requirement comes from 
the fact that the posterior distribution of the polarized in¬ 
tensity P does not follow a normal (Gaussian) distribution. 
Furthermore P is a quantity that must always be positive, 
and this introduces a bias into any estimate. For any true 
Po we would expect to measure on average a polarization 
P > Pq. In order to get the de-biased fluxes, Pdb, from 
the measured ones, P, we choose the Bayesian approach de¬ 


scribed in Vaillancourt (20061 and in Rubino-Martfn et al. 
| |2012b[ ), consisting of integrating the analytical posterior 
probability density function over the parameter space of the 
true polarization. The same posterior can not be used for the 
polarization fraction, 11 = P/I x 100, as it follows a different 
distribution. As, to our knowledge, there is not in the litera¬ 
ture any analytical solution for the posterior distribution of 
n, we numerically evaluate this function by applying Monte- 
Carlo simulations. This approach has already been carried 


| [MT| ). 


out in Lopez-Caraballo et al. (20111 and in Dickinson et al. 


3 RESULTS AND DISCUSSION 
3.1 Maps and consistency tests 

In Fig. we show the intensity map at 11 GHz resulting 
from combining 149 h of observations, where emissions from 
G159.6-18.5, the California nebula (NGC1499) and the 3C84 
quasar are clearly visible. More detailed 7, Q and U maps 
at our four frequencies around the position of G159.6-18.5 
are shown in Fig. The Q and U maps on this source are 
consistent with zero polarization, and therefore upper limits 
on the polarized intensity will be extracted in section |3.3| 
Some stripping is clearly visible in the maps, which is pro¬ 
duced by the presence of regions with a higher noise due to 
a lower integration time per pixel and, to a lesser extent, to 
1// residuals. The inhomogenities in the coverage (integra¬ 
tion time per pixel) maps are caused by the separation of the 
horns in the focal plane, which leads the sky coverage to be 
different when we observe the field before or after crossing 
the local meridian. In the Q and U maps at 11 and 13 GHz of 
Fig-lH the two orthogonal stripes with clearly higher noise 
correspond to regions with integration of ~ 3 — 7 s/pixel 
(pixel size 6.9 arcmin). By comparison, in the central region 
inside the circle where we perform the aperture photometry, 
the integration time is ~ 30 —35 s/pixel, resulting in a lower 
pixel-to-pixel dispersion of the data. 

An important consistency test, that may reveal the 
presence of systematics and other spurious effects, is ob¬ 
tained through jackknife maps. We have uniformly split our 
full dataset in two halves in such a way that the maps of 
number of hits associated to these two halves are as sim¬ 
ilar as possible. The differences of the two halves divided 
by two, for the intensity and polarization maps at our four 
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Figure 2. Intensity and polarization maps at the QUIJOTE four frequency channels around the G159.6-18.5 molecular cloud. The 
Stokes Q and U maps show zero polarization. The solid circle shows the aperture we use for flux integration, whereas the dashed contour 
limits the extent of the ring we use for background subtraction. For comparison, in the last column we plot the corresponding WMAP 
9-year maps at 23 GHz. 


frequencies, are shown in Fig. As expected, the intensity 
emission coming from G159.6-18.5 is consistently cancelled 
out in these maps. While a similar striping pattern to the 
maps in Fig. is seen, these maps are dominated by instru¬ 
mental noise. This is conhrmed by the noise values shown 
in Table where we compare the pixel-to-pixel RMS cal¬ 
culated in two different regions: the external ring that we 
will use for background subtraction when calculating the 
intensity and polarization fluxes, and a region of very low 
sky emission enclosed by the dashed circle represented in 
Fig.H The RMS values in Q and U are similar in the orig¬ 
inal maps and in the sum and difference of the two halves. 
They are typically ~ 250 /rK/pixel (pixel size: 6.9 arcmin) 
or ~ 25 /rK/beam (beam size: 1 ~ degree). In the case of 
the I maps, the RMS in the background ring are higher in 
the sum than in the difference because of the emission of 


the source. In the circle with low sky signal the values in the 
sum and in the difference maps are very similar. 

The last column of Table show the average RMS lev¬ 
els in the Q and U maps normalized by the integration time 
per pixel. As the number of hits per pixel is very inhomo¬ 
geneous, to calculate these numbers we have made a real¬ 
ization of Gaussian noise in which we assign to each pixel 
a noise proportional to , where tpix is the integration 
time per pixel, and then calculate the pixel-to-pixel RMS. 
The amplitudes of the white-noise in the spectra of the time- 
ordered-data range between 898 /rK s^^^ (at 16.7 GHz) and 
1228 /iK s^/^ (at 11.2 GHz;[3 1// residuals make the noises 

^ Note that in this article we are using data from only two out 
of the four horns of QUIJOTE. Then, the global sensitivities of 
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Figure 3. Jackknife difference maps at the four frequencies around the G159.6-18.5 molecular cloud. The structure of this map is 
consistent with instrumental noise, which is predominantly Gaussian (the stripes are associated with regions with a higher noise due to 
a lower integration time per pixel - see text for details). In Tablewe show the pixel-to-pixel RMS values calculated in the background 
annulus enclosed by the two concentric circles of this plot and in the dashed circle. 


calculated on the maps only slightly higher, typically by a 
factor ~ 15%, confirming our previous statement that these 
maps are dominated by white (Gaussian) noise. 

Another important consistency test for the presence of 
systematics, and in particular for the / to Q/U polariza¬ 
tion leakage, is to confirm that our polarization maps are 
consistent with noise in the position of unpolarized sources. 
This verification is provided by the nearby California HIT 
region, which is also covered by our observations. As any 
standard HIT region, it is dominated by free-free emission at 

the experiment are a factor better, i.e. ^ 650 /rK the 

number quoted in Section [2Tj 


the QUIJOTE frequencies, which is known to be practically 
unpolarized. In Fig. we show intensity and polarization 
maps towards the California region, showing that the Q and 
U maps are consistent with zero polarization. 


3.2 Intensity SED 

As was indicated in section |2.2| we take the COSMOSO- 
MAS fluxes for G159.6-18.5 from IPlanck Collaboration et] 
E] ( |2011[ ), which have already been corrected for the flux 
loss caused by the filtering of COSMOSOMAS data. While 
here we will use aperture photometry to derive our fluxes, in 
[Planck Collaboration et al.| (|2014a|) they were obtained by 
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Figure 4. Intensity and polarization maps around the California HII region, at 17 GHz from QUIJOTE and at 23 GHz from WMAP 
9-year data. The Stokes Q and U maps show zero polarization as is expected for a free-free dominated region. This test therefore 
demonstrates the absence of significant polarization systematics inherent to our experiment or to the data processing. The ellipses show 
the aperture and the background annulus that we use to derive polarization upper limits in this region. 


u (GHz) 


(/iK/pixel) 

o'Q 

(/iK/pixel) 


(/iK/pixel) 

CQ.U (mK 8^/2) 


Map 

Sum 

Dili. 

Map 

Sum 

Diff. 

Map 

Sum 

Diff. 

Map 

11.2 

1592 

1593 

1277 

389 

396 

394 

361 

361 

361 

1478 

12.9 

1328 

1377 

1119 

300 

308 

314 

293 

301 

297 

1179 

16.7 

861 

861 

755 

167 

169 

171 

166 

166 

161 

1158 

18.7 

1040 

1041 

966 

199 

202 

213 

201 

205 

201 

1461 

11.2 

1377 

1415 

1320 

406 

418 

460 

315 

335 

346 

1283 

12.9 

1331 

1512 

1529 

352 

387 

349 

309 

339 

351 

1192 

16.7 

672 

673 

651 

142 

142 

140 

125 

125 

136 

1009 

18.7 

803 

809 

876 

149 

150 

157 

164 

165 

163 

1413 


Table 1. RMS pixel-to-pixel (pixel size 6.9 arcmin) calculated in the original maps, and in the average and difference divided by two of 
the two jackknife maps, in two different regions. The values above and below the horizontal line have been calculated respectively in the 
background annulus and in the dashed circle depicted in Fig.[^ The numbers in the last column correspond to the average noises in the 
original Q and U maps normalized by the integration time per pixel (units: fiK s^/^). 
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Frequency 

(GHz) 

Flux 

density (Jy) 

Flux density 
residual (Jy) 

Telescope/ 

survey 

0.408 

10.5 ±4.0 

2.7 ± 4.2 

Haslam 

0.82 

7.4 ± 2.1 

0.1 ± 2.4 

Dwingeloo 

1.42 

6.8 ± 1.6 

-0.1 ± 2.0 

Reich 

10.9 

16.1 ± 1.8 

10.3 ± 1.9 

COSMO. 

11.2 

15.0 ±2.0 

9.2 ± 2.1 

QUIJOTE 

12.7 

20.0 ±2.2 

14.2 ±2.2 

COSMO. 

12.9 

18.1 ± 2.1 

12.4 ± 2.1 

QUIJOTE 

14.7 

28.4 ±3.1 

22.6 ±3.1 

COSMO. 

16.3 

35.8 ±4.0 

30.0 ±4.0 

COSMO. 

16.7 

33.9 ±2.4 

28.1 ± 2.5 

QUIJOTE 

18.7 

35.2 ±3.7 

29.2 ±3.8 

QUIJOTE 

22.8 

40.2 ±2.4 

34.0 ±2.5 

WMAP 

28.4 

40.4 ±2.4 

33.5 ±2.6 

Planck 

33.0 

38.1 ± 2.4 

30.4 ± 2.8 

WMAP 

40.7 

32.8 ±2.5 

23.1 ± 3.3 

WMAP 

44.1 

29.8 ±2.6 

19.1 ± 3.7 

Planck 

60.8 

27.5 ±3.8 

8.3 ±6.9 

WMAP 

70.4 

32.3 ±4.9 

5.2 ± 10.0 

Planck 

93.5 

59.5 ±9.3 

1.8 ± 22.3 

WMAP 

100 

81 ± 17 

11 ± 10 

Planck 

143 

194 ± 24 

-17± 82 

Planck 

217 

1011 ± 122 

196 ± 320 

Planck 

353 

4286 ± 446 

344 ± 1376 

Planck 

545 

14858 ± 1470 

208 ± 4588 

Planck 

857 

45235 ± 4045 

-1352 ± 13168 

Planck 

1249 

86696 ± 6674 

-4878 ± 25315 

DIRBE 

2141 

114650 ±6891 

6845 ± 43590 

DIRBE 

2998 

54361 ± 2624 

-837 ± 28264 

DIRBE 


Table 2. Flux densities for G159.6-18.5 in the Perseus molecular 
cloud. All fluxes have been calculated through aperture photom¬ 
etry in a ring of radius 1.7° and subtracting the median of the 
background in a ring between 1.7° and 1.7\/2 , except those com¬ 
ing from the COSMOSOMAS experiment (10.9, 12.7, 14.7 and 
16.3 GHz), which have been taken from [Wanck Collaboration et| 
|al.| ) |20lT| l. Also shown are the residual AME fluxes, obtained after 
subtraction of the free-free, CMB and thermal dust components. 
The last column indicate the telescope or survey from which the 
data have been extracted. 


fitting the amplitude of an elliptical Gaussian with a fixed 
size of 1.6° X 1.0° (FWHM). In a first-order approximation, 
the fluxes obtained through Gaussian fitting will be equiv¬ 
alent to those obtained from aperture photometry using a 
given aperture size. Therefore, in order to get a reliable in¬ 
tensity SED we choose a size for the aperture that gives the 
most similar fluxes to those presented in [Planck Gollabora- 


|tion et al.| (|201l1 for the|Haslam et al.| | |1982 |, Berkhuijsen 
1 1972[ |, [Reich fc Reich (19861, WMAP , Planck and DIRBE 
maps (it must be noted that the WMAP and Planck maps 
used in this work correspond to a different release to that 
used in |Planck Collaboration et al. (20111, but this will have 
a negligible effect). After trying different values, we found 
that a radius of 1.7° gives the best agreement, with a very 



Figure 5. Spectral energy distribution of G159.6-18.5. QUIJOTE 
data points are depicted in red, together with other ancillary data 
including COSMOSOMAS, WMAP 9-year data and Planck data. 
At intermediate frequencies the excess emission associated with 
the AME clearly shows up. A joint fit has been performed to all 
the data points except 100 and 217 GHz (which are affected by 
CO contamination), consisting in the following components; free- 
free (dotted line), spinning dust from a mixed environment with 
a high-density molecular (long-dashed line) and a low-density 
atomic (dashed line) gas, CMB (dash-dotted line) and thermal 
dust (dashed-triple-dotted line). The solid line represents the sum 
of all the components. 


low reduced chi-squared of Xred ~ 0.098. The median of the 
background is computed in an external ring between 1.7° 
and 1.7\/2 , which has the same area as the aperture. The 
derived fluxes in the QUIJOTE maps and in the other an¬ 
cillary maps are listed in Table 

The final SED is depicted in Fig.[^ where the presence 
of AME clearly shows up at intermediate frequencies as an 
excess of emission over the other components. The inten¬ 
sities derived from these QUIJOTE observations trace, for 
the hrst time after the original measurements of the COS¬ 
MOSOMAS experiment (Watson et al. 20051, the down¬ 
turn of the spectrum at frequencies below ~ 20 GHz, as 
it is predicted by spinning dust models. In total, 13 data 
points are dominated by AME: the four QUIJOTE points, 
the four COSMOSOMAS points, the WMAP 22.8, 33.0 and 
40.7 GHz frequencies and the Planck 28.4 and 44.1 GHz fre¬ 
quencies. We perform a joint multi-component fit to all the 
data, consisting of free-free emission, which dominates in 
the low-frequency tail, spinning dust, which dominates the 
intermediate frequencies, a CMB component, and thermal 
dust, which dominates the high-frequency end. As it was 
done in Planck Collaboration et al. (20111, to avoid possible 
CO residuals we exclude from the fit the 100 GHz and the 
217 GHz values. We fix the spectrum of the free-free using 
the standard formulae shown in IPlanck Collaboration et al.l 
(20111, with a value for the electron temperature typical of 
the solar neighbourhood, Te = 8000 K, and fit for its ampli- 
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tude, which is parameterized through the emission measure 
{EM). _ 


Following Planck Collaboration et al. ( 20111 , we con¬ 
sider a high-density molecular phase and a low-density 
atomic phase which produce spinning dust emission, and ht 
their respective amplitudes, which are given by the hydro¬ 
gen column densities a-nd A^h*- The CMB amplitude 

is denoted by ATcmb, and would correspond to the average 
of the primordial CMB fluctuations in the aperture. Finally, 
the thermal dust is modelled as a single-component modi- 
hed blackbody curve, T 25 o{v/ 12 QQ GYiz)^'^Bi,{Td), which de¬ 
pends on three parameters: the optical depth at 250 /rm 
(''" 250 ), the emissivity spectral index (/3d) and the dust tem¬ 
perature (Td). Therefore, we jointly fit 7 parameters to all 
the data points: EM, N^°\ ATcmb, T250, / 3 d and Td. 

To define the spinning dust spectra of t he molecular and 
atomic phases we use the SPDUST.2 cod^ (Ali-Haimoud et 


aL]|2009 Silsbee et al.||2011|. Initially we use the same val 


Planck Collaboration et al. (20111 for the different 


parameters on which the spinning dust emissivity depends. 
This involves a modification of the code, as by default it 
uses ao,i = 0.35 nm and ao ,2 = 3.0 nm for the centroids 
of the two lognormal functions defining the dust grains size 
distributions, while the carbon abundance per hydrogen nu¬ 
cleus, be, is selected from any of the values in table 1 of 


Weingartner & Draine (20011. Following Planck Collabo¬ 


ration et al. (20111, we use instead ao.i = 0.58 nm and 


0 - 0,1 = 0.53 nm, respectively for the molecular and atomic 
phases, and be = 68 ppm. Using these two models, we get a 
good fit to our full dataset, with Xred ~ 0.99, where the spin¬ 
ning dust component is clearly dominated by the molecular 


phase, as was found by Planck Collaboration et al. (20111. 


To analyze the possibility of the existence of slightly 
different spinning dust models that could provide a better 
ht to the data, we take the AME residual fluxes from the 
previous ht, and produce a grid of models varying some of 
the parameters of the molecular phase component. In par¬ 
ticular, we vary: i) the hydrogen number density tih between 
10 and 500 cm“^ with a step of 5 cm“^; ii) the kinetic gas 
temperature Tg between 5 and 200 K in steps of 5 K; iii) 
the intensity of the radiation held relative to the average 
interstellar radiation held, for which we consider only the 
values Go = 1 and 2; iv) and the hydrogen ionization frac¬ 
tion, for which we consider the values xh = 10, 112, 1000 
and 10000 ppm. The best ht is obtained for Go = 1 and 
Xu = 112 ppm, the same values used in [Planck CoIlabora-| 
|tion et ^ ( |201l| . As the ht is very degenerate, to dehne the 
most-likely values for nu and Tg we set Gaussian priors on 
four diherent parameters. First, we put soft priors on nu and 
Tg cen tred on the same values used in [Planck Collaborationj 


et al. 


(20111, (nujo = 250 cm“^ and (Tg)o = 40 K, and with 
standard deviations (j{nu) = 80 cm“® and cr(Tg) = 60 K. 
An additional prior on can be derived from the canon- 
H cm-2 


ical relation 2.13 x 10^ 


= Ixioo (Finkbeiner et 


ar]|2004 1. Using T 2 so and /3d from our best-ht model of the 
thermal dust component, we extrapolate the optical depth 
to 100 /rm, and hnd Nu°'^ = 2.909 x lO^'^ H cm“^. We use 
this value to dehne the centre of the Gaussian prior, and 
= 2 X 10^^ H cm“^. Finally, the ratio between the 
hydrogen column density and the hydrogen volume density, 
2 = /nu, gives an estimate of the length along the line 
of sight of the spinning-dust-emitting region. We assume 
that this length might be of the order of the transverse an¬ 
gular size of the source. The source subtends an angle of 
around 2°, which at the distance of the Perseus complex, 
260 pc (Cernicharo et al. 19851, corresponds to 9.08 pc. 


We therefore set a fourth prior on this quantity defined by 
zq = 9.08 pc and a{z) = 4 pc. 

In Fig. we show the marginalized likelihoods over nu 
and Tg. We define the best-values for these parameters from 
the 50% integrals of the probability distribution, and the 
confidence intervals from the region encompassing the 68% 
of the area around those values. We get nu = 223.2/'/g2 g 
and Tg = 60.2'/t^2 g K. As it was said before, the values of 
the intensity of the radiation field and of the hydrogen ion¬ 
ization fraction that maximize the likelihood are Go = 1 
and Xu = 112 ppm, respectively. We fix the other param¬ 
eters of the molecular phase, and all the parameters corre¬ 
sponding to the atomic phase, to the same values that were 


used in Planck Collaboration et al. (20111. All these values 
are shown in Table In this table, xc represents the ion¬ 
ized carbon fractional abundance, y the molecular hydrogen 
fractional abundance, and j3 the average dipole moment per 
atom. The meaning of the other parameters have been ex¬ 
plained before in the text. We then obtain the corresponding 
spinning dust spectra for the molecular and atomic phases 
using these parameters as inputs for SPDUST.2. Fixing these 
spectra, we perform a joint fit of the hve aforementioned 
components, obtaining the best-fit values for the 7 parame¬ 
ters defining these models, which are also shown in Table 
We get Xred = 0.99, the same value as before, so we do not 
manage to improve the quality of the global fit after improv¬ 
ing the spinning dust models. This highlights the difficulty 
of constraining the parameters on which the spinning dust 
emission depends due to the the strong degeneracies between 
them. 

The total hydrogen column density is (4.40 ± 0.83) x 
lO^'^ H cm“^. This is a bit higher than the expected value 
of 2.89 X 10^^ H cm~^, which has been derived from the 
aforementioned tioo — -^h canonical relation, and extrapo¬ 
lating Tioo to r 25 o using the fitted spectrum for the thermal 
dust emission. The inferred lengths of the two spinning dust 
phases along the line of sight, 2 = Nu/nu, are respectively 
2 ™°' = 5.63T0.76 pc and 2 ^^* < 12.6 pc (a 68.3% C.L. upper 
bound is used here, as the error bar is higher than the es¬ 
timate). These values are of the order, or compatible, with 
the transverse size of the region, 9.08 pc. Our htted values 
for /3d and Td are consistent within 1-sigma with those de¬ 


rived in Planck Collaboration et al. (20111. On the other 


http: / / www.tapir.caltech.edu/~yacine / spdust / spdust.html 


hand, we get lower values for EM, Nu and r 25 o, but this 
is because these depend on the solid angle subtended by 
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Figure 6. Likelihood distributions for the hydrogen number density (left) and for the kinetic gas temperature (right) associated with 
the molecular phase of the spinning dust component. These were obtained after marginalizing over the rest of parameters on which the 
spinning dust emissivity depends (see text for details). The vertical lines show the most likely value, defined through the 50% integral of 
the cumulative probability distribution, and the dashed regions depict the 68% confidence interval around this value. 


EM (cm ®pc) 

23.9 ± 1.9 


Molecular 

Atomic 

nil (cm“®) 

223.2 

30 

Go 

1 

2 

Tg (K) 

60.2 

100 

(ppm) 

112 

410 

xc (ppm) 

1 

100 

y 

1 

0.1 

ao,i (nm) 

0.58 

0.53 

ao ,2 (nm) 

3.0 

3.0 

be (ppm) 

68 

68 

/3 (D) 

9.34 

9.34 

Ah (lO^i cm-2) 

3.88 ±0.52 

0.52 ±0.65 

ATcmb (yK) 

22.6 ± 13.6 

dd 

1.73 ± 0.11 

Td (K) 

18.2 ± 0.6 

T250 

(2.78 ± 0.41) X 10-4 

X?ed 

0.99 


Table 3. Model parameters. The emission measure EM is fitted 
to the data and defines the amplitude of the free-free emission. 
Parameters from nn to /3 (see section |3.2| for an explanation of 
their meaning) are used as inputs for SPDUST.2, and define the 
shape of the spinning dust spectra for the molecular and atomic 
phases. The amplitude of these spectra are determined through 
the fit to the data, and are given by the hydrogen column density 
Nil. The best-fit values for the amplitude of the CMB component 
ATcnib for the three parameters defining the thermal dust 
spectrum, as well as the reduced chi-squared, are also shown. 


the region. This value is a factor « 5 higher in our case, 
because we are using aperture photometry instead of Gaus¬ 
sian fitting. When this factor is taken into account, then 
our values are brought into a better agreement with those 


of Planck Collaboration et al. (2011 b Finally, we find a post 


five value for ATcmb, whereas in |Planck Collaboration et al.| 

| |2011[ | it was negative. However, we have checked that our 
value agrees with the average level of the CMB anisotropies 
within the aperture, which has been found to be 23.6 fiK 
in the Planck-DKl CMB map resulting from the SMICA 


component separation methocj^ (|Planck Collaboration et al. 
2014d|. 


In Fig. we show the residual spinning dust spectrum, 
obtained after subtracting the best-fit free-free, CMB and 
thermal dust components. 


3.3 Polarization constraints 


As no clear emission is seen in the Q and U maps of Fig.[^ we 
derive upper limits on the polarization fraction of G159.6- 
18.5, following the procedure explained in section [231 In sec¬ 
tion we defined the sizes of the aperture and of the back¬ 
ground annulus so that we reproduced the fluxes in |Planck| 
Collaboration et al. ( 2011| , which were obtained through 
Gaussian fitting. In order to minimize the error associated 
with the background subtraction, here we extend the size of 
the background annulus, and use a circular aperture with ra¬ 
dius 1.5°, and a background ring between 1.5° and 2.5°, with 
an extension towards the north-west as is shown in the maps 
of Fig. 1^ At each frequency we calculate the RMS levels in 
this background annulus, to define the quantities o-{Tj) that 
are introduced in equation to calculate the errors on the 


^ Downloaded from the Planck Legacy Archive, 
http: //www. sciops. esa. int/index .php?project=planck.!cpage= 
Planck_Legacy_Archive. 
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Figure 7. Spectrum of the AME in G159.6-18.5 after subtract¬ 
ing the best-fit free-free, CMB and thermal dust components. 
The long-dashed line and the dashed line show respectively the 
spinning dust spectrum for the molecular and the atomic phases, 
while the solid line is the sum of the two. 


Stokes parameters Q and U. These errors, together with the 
fluxes resulting from the aperture photometry integration, 
are quoted in Table Note that here the intensity fluxes 
in the QUIJOTE frequencies are slightly different to those 
presented in Table [Rowing to the different apertures. In or¬ 
der to get the AME residual fluxes shown in Table using 
the new fluxes, we repeat the same fit that was performed 
in section |3.2| considering the same spinning dust spectra 
for the molecular and atomic phases. 

It is important to point out that, as it became clear 
when we discussed the results of the jackknife tests in sec- 
tion |3.1[ the uncertainties on the I fluxes are here biased high 
because some extended emission of the source leaks into the 
background annulus. This does not have significant implica¬ 
tions in our analysis as the uncertainty in the polarization 
fraction is driven by the errors in Q and U. 

The Q and U fluxes shown in Table are consistent 


with zero, and therefore we obtain de-biased (in section 2.3 


we explained how this de-biasing is applied) upper limits 
at the 95% confidence level on the polarized intensity, Pdb- 
We also show in Table upper limits at the 95% confidence 
level on the polarization fraction, taking as reference both 
the total intensity (Ildb) and the residual AME intensity 
(IlAME.db). The values below the horizontal line in this table 
correspond to constraints obtained in maps that have been 
built by combining the two frequency bands of each horn. 
The most stringent upper limit we get on the polarization 
fraction is Ildb < 2.85%, and is obtained after combining 
the maps at 16.7 and 18.7 GHz. Due to the decrease of the 
intensity flux at lower frequencies, the constraints at 11.2 
and 12.9 GHz are less stringent. 

Note that, under the reliable assumption that the free- 
free emission is unpolarized | Rybicki fc Lightman|1979 1, any 
possible detection of polarization below i/ ~ 30 GHz, where 
the thermal dust is clearly sub-dominant, should in principle 
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be ascribed to AME. One caveat to this hypothesis is the 
possible presence of a Faraday Screen (FS) hosting a strong 
and regular magnetic field, which could induce a rotation 
of the background polarized emission. This idea was pro¬ 


posed by Reich & Reich (20091 to explain the high degree of 


polarization they detected towards G159.6-18.5 in observa¬ 
tions from the Effelsberg telescope at 2.7 GHz. They suggest 
that this same mechanism could indeed be the responsible 
for the tentative polarized emission seen by |BattisteIli et ah] 
(20061 at 11 G Hz. Using as face valu e the rotation mea¬ 
sure obtained by Reich & Reich (20091, RM= 190 rad m“^. 


Lopez-Garaballo et al. (20111 estimated a polarization frac¬ 


tion of < 0.2%, well below the upper limit at this frequency. 
Using this RM we estimate polarization fractions from the 
FS of « 3.5% and « 0.4%, respectively at 11 and 19 GHz. 
These values are well below our upper limits, while the value 
at 11 GHz is compatible with the measurement of |Battistelli| 
et al. (20061. 


In section [3T] we mentioned that one important consis¬ 
tency test for our data processing is the verification that it 
shows no polarization at the position of un-polarized sources 
like the California HII region. In Table we show the corre¬ 
sponding I, Q and U fluxes, and derived upper limits on P 
and n, which have been obtained using the apertures shown 
in Fig. The limits on the polarization fraction stand at 
the level of 2 to 4% (95% C.L.), depending on the frequency 
band. This constrains the possible existence of a leakage 
from intensity to polarization to values below this level. In 
fact, as it was mentioned in section [Hi we expect the leak¬ 
age to be below 0.3% (^ —25 dB). This number has been 
verified using observations of Gas A, where we recover a 
leakage pattern with a quadrupolar shape in the Q and U 
maps, with a level of around 0.2% (further details will be 
presented in separate technical papers). 

We plot our constraints, and previous results in the lit¬ 
erature at different frequencies, in Fig. It can be seen 
here that our observations fill the gap between previous re¬ 
sults at frequencies below 11 GHz and above 20 GHz. The 
point at 9.65 GHz is an upper limit coming from GBT 


observations on LDN1622 (Mason et al. 20091. The value 


at 11 GHz represents the aforementioned tentative detec¬ 
tion towards G159.6-18.5 at the 1.8o' level (H = 3.4)^l g%) 
of AME polarization from the COSMOSOMAS experiment 


(Battistelli et al. 20061, which is still consistent with our 


upper bound at 12.0 GHz. The points at iz > 20 GHz come 


from WMAP 7-year data on the Perseus (Lopez-Caraballo 


et al. 20111 and p Ophiuchi (Dickinson et al. 2011, they 
also present results on Perseus) molecular clouds and on 
the HII region [LPH96] 201.663-1-1.643 (Rubino-Martln et 
al. 2012b, they also show less-stringent constraints from 
the Pleiades reflection nebula and from the dark nebula 
LDN1622). Currently, the most stringent constraint is that 
obtained by [Lopez-Caraballo et ab ( 2011| ): H < 0.98% at 
23 GHz. This limit benefits from the fact that it is obtained 
at a frequency close to the peak (« 28 GHz) of the AME 
in G159.6-18.5. The constraints on the polarization fraction 
from QUIJOTE come from a spectral region where the in- 
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iz (GHz) 

I (Jy) 

7ame (Jy) 

Q (Jy) 

U (Jy) 

Pdb (Jy) 

Hdb (%) 

llAME.db (%) 

11.2 

14.0 

± 

3.3 

9.4 

± 

3.4 

-0.05 

± 

0.59 

-0.39 

± 

0.47 

< 

1.19 

< 

8.79 

< 13.85 

12.9 

17.5 

± 

3.5 

12.5 

± 

3.5 

0.45 

± 

0.68 

-0.23 

± 

0.56 

< 

1.44 

< 

8.46 

< 12.20 

16.7 

31.2 

± 

3.0 

28.2 

± 

3.0 

0.24 

± 

0.47 

-0.83 

± 

0.52 

< 

0.95 

< 

5.14 

< 5.71 

18.7 

31.5 

± 

4.6 

29.4 

± 

4.7 

-0.11 

± 

0.75 

0.18 

± 

0.70 

< 

1.49 

< 

4.76 

< 5.08 

12.0 

15.5 

± 

3.4 

10.2 

± 

3.4 

0.14 

± 

0.38 

-0.36 

± 

0.37 

< 

0.94 

< 

6.26 

< 10.10 

17.7 

31.8 

± 

3.7 

26.7 

± 

3.7 

0.06 

± 

0.42 

-0.26 

± 

0.40 

< 

0.91 

< 

2.85 

< 3.42 


Table 4. Intensity and polarization fluxes at the position of G159.6-18.5 determined through an aperture photometry integration, and 
upper limits on the polarized intensity and on the polarization fraction. We show the total intensity fluxes (/) as well as the residual 
AME fluxes (/ame) after subtracting the free-free and thermal dust components. Similarly, for the polarization fraction we show upper 
limits calculated with respect to the total and to the residual intensities. The error bars are at 1-sigma, whereas the upper limits are 
given at the 95% confidence level. Results are shown for each of the four QUIJOTE frequency channels, as well as for a combination of 
the low and high frequency bands. 


u (GHz) 

I (Jy) 

Q (Jy) 

U (Jy) 

Pdb (Jy) 

Hdb (%) 

11.2 

59.0 ± 4.1 

-0.45 ±0.36 

0.52 ± 0.42 

< 1.25 

< 2.12 

12.9 

54.7 ± 4.2 

-0.02 ±0.44 

0.24 ± 0.46 

< 0.96 

< 1.76 

16.7 

49.0 ±3.7 

-0.07 ± 0.38 

0.38 ±0.35 

< 0.91 

< 1.87 

18.7 

51.8 ±6.2 

-0.10 ± 0.52 

-0.68 ± 0.50 

< 1.43 

< 2.82 


Table 5. Intensity and polarization fluxes of the California HII region determined through an aperture photometry integration, and 
upper limits on the polarized intensity and on the polarization fraction. The error bars of the I, Q and U fluxes are at 1-sigma, whereas 
the upper limits on the polarized intensity and on the polarization fraction are referred to the 95% confidence level. 


tensity flux density drops. However, according to the model 


of Hoang et al. (20131, the polarization fraction of the spin¬ 
ning dust emission peaks at « 5 GHz, and therefore observa¬ 
tions at lower frequencies are more appropriate to constrain 
this model. All measurements, including ours, are consis¬ 
tent with this model, except maybe the previous limit from 
|L6pez-Caraballo et al.| ( |2011[ ) and that from [Mason et al.| 
(20091, which are slightly below. 


Our upper limits on H are not only affected by the de¬ 
crease of AME flux at our frequencies, but also by the sparse 
sky coverage of our observations. We decided to survey an 
area large enough to cover the California HII region and also 
to get the source simultaneously in the four horns. Currently 
we are undertaking observations on a smaller area, centred 
in two individual horns, with the goal to increase by a fac¬ 
tor 8 the integration time per unit area. This would improve 
our map sensitivity by a factor 2.8, which would help push 
our current upper limits down to ~ 1.7% at 18.7 GHz, and 
« 1.0% after combining the 16.7 and 18.7 GHz frequency 
bands, which is well below the polarization degree predicted 


by the Hoang et al. (20131 model at these frequencies. Fur- 


TGI experiment, which is expected to be 13 times more 
sensitive than the current MFI, will push these upper limits 
down by a significant amount. 

However, it is important to note that the curve of the 


polarization degree obtained by Hoang et al. (20131 rep¬ 


resents an upper limit. They inferred the alignment effi¬ 
ciency of interstellar dust grains from observations of the 
UV polarization excess towards two stars to predict the 
polarization degree of the spinning dust emission. In so 
doing they assumed that the UV polarization bump is 
produced exclusively by polycyclic aromatic hydrocarbon 
(PAHs) molecules, those which are thought to be respon¬ 
sible for spinning dust emission. However, if the graphite 
grains were also aligned, contributing to the UV polarization 
bump, then the alignment efficiency of PAHs could actually 
be lower and so would be the inferred degree of spinning 
dust polarization. 


thermore, future observations at 30 GHz with the upcoming 


The previous model from Lazarian & Draine (20001 for 
spinning dust polarization is also plotted in Fig. for the 
case of a cold neutral medium (CNM), together with dif¬ 
ferent models from Draine & Lazarian ( |1999 1 predicting the 
degree of polarization of the magnetic dipole emission for dif- 
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Frequency (GHz) 

Figure 8. Constraints (95% C.L.) on several microwave emission 
models based on our upper limits on G159.6-18.5 (black dots) 
and on others in different regi ons from the literature ||Mason et 
al. 12009] [Battistelli et al.|2006||L6pez-Caraballo et al.|2()ll||Dick- 

inson et al.||2011| jRubino-Martm et al.||2012b| |. The black lines 
are theoretical predictions for the frequency dependence of the 
fractional polarization of the electric dipole (ED) emission from 
jLazarian &: Drainej | |2000| | (dashed line) and jHoang et ar] | |2013[ l 
(solid line). The red and blue lines correspond to the models of 
jDraine fc Lazarian| ( |199'9t for the magnetic dipole (MD) emission. 
The different lines correspond to different grain compositions and 
geometries. The observations are generally consistent with the 
electric dipole emission models, and rule out the models for the 
fractional polarization of the magnetic dipole emission. 


ferent grain shapes and compositions. In particular, we show 
the cases of grains made of metallic Fe, and of the hypothet¬ 
ical material X4 defined in [Draine fc Lazarian] ( |1999| ) . Con¬ 
cerning the grain geometries, different axial ratios al:a2:a3 
are shown, as indicated in the figure. All these models pre¬ 
dict polarization fractions typically above 10%, much higher 
than that of the spinning dust emission, and are ruled out 
by our and previous observations. However, it must be noted 
that these models consider the dust grains to have a perfect 
ordering of the atomic magnetic moments in a single do¬ 
main. A more realistic case considering randomly-oriented 
single-domain magnetic inclusions has been recently studied 
by Draine & Hensley (2013), resulting in much lower polar¬ 
ization degrees, of the order of 5% in the range 10 — 20 GHz, 
or below. 


4 CONCLUSIONS 

We have presented the first results obtained with the QUI¬ 
JOTE experiment, a new polarimeter aimed at measuring 
the B-mode anisotropy from inflation and also at charac¬ 
terizing the polarized foregrounds at low frequencies. These 
results are based on observations of the Perseus molecular 
complex, one of the regions where the AME has extensively 
been studied. Our observations cover G159.6-18.5, the region 


where AME is produced, and the Galifornia HII region. Our 
intensity data on G159.6-18.5 traces the decrease of the flux 
of this source at frequencies below ~ 20 GHz, confirming the 
prediction of the models based on spinning dust radiation, 
which currently are the best at reproducing the frequency 
spectrum of the AME. This confirms previous results on this 


region obtained with the COSMOSOMAS experiment | Wat- 
son et aL]|2005 |. When combining the QUIJOTE measure¬ 


ments with data coming from GOSMOSOMAS, WMAP-9yr 
and Planck-DKl, we get an intensity SED with 13 points be¬ 
tween 10 and 50 GHz being dominated by AME, providing 
what probably is the most precise spectrum of this emission 
ever measured in an individual object. This intensity SED is 
well fitted (Xred = 0.99) by a combination of free-free, CMB, 
spinning dust and thermal dust, and two spinning dust com¬ 
ponents associated to a high-density molecular phase and to 
a low-density atomic phase. We attempt to fit some of the 
parameters describing the physical environment of the cloud, 
and which define the shape of the spinning dust spectrum, 
but the solution is clearly degenerate and this approach does 
not achieve a better fit. However, with the physical param¬ 
eters we use for the spinning dust spectra, we get reliable 
values for the 7 parameters of our model. 

The California HII region is well detected in intensity, 
presenting a spectrum at QUIJOTE frequencies consistent 
with free-free emission. The polarization maps at the posi¬ 
tion of California are consistent with noise, as it is expected 
for regions dominated by free-free emission, which is very 
lowly polarized. No polarization is detected at the position 
of G159.6-18.5 either, so then we extract upper limits on 
the polarized intensity and on its degree of polarization. Af¬ 
ter combining the low-frequency and high-frequency maps of 
two horns of QUIJOTE, we get H < 6.26% and H < 2.85% 
(95% C.L.), respectively at 12.0 and 17.7 GHz, being the 
first upper limits on the AME polarization fraction in this 
frequency range. These constraints are consistent with mod¬ 
els of the spinning dust emission, that predict polarization 
fractions between 2 and 3% in this frequency range, and 
rule out several models based on magnetic dipole emission. 
Pushing these constraints even further is crucial to test the 
theoretical models and also to assess up to what level AME 
will hinder the detection of the B-mode signal by current 
and future experiments operating in the frequency range 10 
to 80 GHz. New QUIJOTE observations in this region will 
improve the sensitivity in our maps so that we will be able 
to improve the previous constraints at least by a factor 2, 
allowing to reach a sensitivity on the polarization fraction 
close to the predictions of the spinning dust models. 
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